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Abstract. We study the relation between the measured anisotropies in the extragalactic 
diffuse gamma-ray background (DGRB) and the DGRB spectral intensity, and their potential 
origin from the unresolved blazar population. Using a physical-evolution model for blazars 
with a luminosity dependent density evolution (LDDE) and an observationally-determined 
luminosity-dependent blazar spectral energy distribution (SED), we find that blazars can ac- 
count for the observed anisotropy of the DGRB consistent with their observed source-count 
distribution, but are in turn constrained in contributing significantly to the observed DGRB 
intensity. For the best-fit LDDE model accounting for the DGRB anisotropy and source- 
count distribution, blazars only contribute 5.7^'q% (68% CL) of the DGRB intensity above 
1 GeV. Requiring a higher fraction of the DGRB intensity contribution by blazars over- 
produces the DGRB anisotropy, and therefore blazars in the LDDE+SED-sequence model 
cannot simultaneously account for the DGRB intensity as well as anisotropy. We discuss 
the limitations of LDDE models. However, these models do not require the many unjustified 
and observationally-inconsistent simplifying assumptions — including a single power law for all 
blazar spectra and a simple broken power-law model for their source-count distribution — that 
are present in much previous work. 
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1 Introduction 

The launch of the Large Area Telescope (LAT) aboard the Fermi Gamma Ray Space Tele- 
scope [1] has ushered in a new era of an accurate understanding of the extragalactic gamma- 
ray sky. Blazars are the most numerous extragalactic gamma-ray sources seen by Fermi-LAT, 
yet the contribution of blazars to the extragalactic gamma-ray sky remains an open question. 
Below the Fermi-LAT point-source resolution threshold, the inferred total gamma-ray flux in 
the diffuse component due to blazars is highly model dependent, and is also every dependent 
on the nature of the blazar population with the spectrally-dependent ability of Fermi-LAT 
to resolve point sources. 

Blazars are comprised of two classes of active galactic nuclei ( AGN) : flat-spectrum radio 
quasars (FSRQs) and BL Lacertae objects (BL Lacs). Both of these sources have highly 
variable luminosity and high luminosity [1], with a bolometric luminosity which is dominated 
by their gamma-ray luminosity [2-4], so they are believed to come from a single source class, 
blazars. A blazar occurs when an AGN is observed down the relativistic jet rather than the 
usual observation of the accretion disk surrounding the AGN. Unlike the near-isotropic flux 
coming from the accretion disk, the emission from the jet is expected to be relativistically 
beamed [5, 6]. 

Several models of the blazar population have been proposed in the literature [7-22]. 
Physical-evolution models are formed by considering the blazar evolution and their con- 
stituent spectral dependence. A pure luminosity evolution model of blazars considers only 
the evolution of blazar luminosity with redshift [10, 12, 16, 18]. Luminosity-dependent den- 
sity evolution (LDDE) models also allow blazar peak luminosity to vary with redshift [23, 24]. 
These models are based on the observed properties of X-ray AGN and relate the gamma-ray 
properties of the blazar model to these X-ray AGN measurements. Some blazar models, such 
as those from Ref. [19] (hereafter FB10) and Ref. [25] (hereafter CKS), assume that the mea- 
sured Fermi-LAT blazar spectral index distribution is largely independent of blazar flux and 
that the measured Fermi-LAT source-count distribution is a power law that continues unal- 
tered down to zero flux. In addition, the spectral properties of blazars are often assumed to 
be a simple power-law or a distribution of power-laws, but more detailed frequency-dependent 
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models have been used as well [14, 20, 26]. These frequency-dependent blazar spectra provide 
a more realistic fit to measured blazar properties, seen in the spectral energy distribution 
(SED) sequence of Ref. [2-4]. 

In this work, we study the LDDE with a SED sequence blazar model like that in Ref. [20] 
(hereafter ABH). The use of an LDDE blazar evolution model incorporates explicitly the 
blazar evolution in both redshift and luminosity. The blazar SED sequence gives the energy- 
dependence of the blazar luminosity and accounts for differing spectra for blazars of different 
luminosity. Only such a physically-motivated model of the blazar evolution and spectrum 
can accurately predict the blazar flux below the Fermi-LAT threshold. In fact, the measured 
extragalactic diffuse gamma-ray background (DGRB) by the Fermi-LAT Collaboration [27] 
was predicted in shape and amplitude by Inoue & Totani [26] using the LDDE plus SED 
sequence blazar evolution model [28]. 1 

In particular, the properties of the blazar evolution and spectrum strongly affect the 
blazar contribution to the DGRB. The DGRB is derived by measuring the full gamma-ray 
spectrum, then removing resolved point-sources and a model for the Galactic diffuse flux. The 
remaining background, the DGRB, is made up of unresolved extragalactic point-sources, such 
as blazars, and the isotropic component of unmodeled Galactic sources, such as millisecond 
pulsars or, potentially, annihilating dark matter. The most accurate measurements of the 
DGRB have been made with the Fermi-LAT [19, 27, 29, 30]. 

The DGRB spectrum measured by the Fermi-LAT has an intensity /(> 100 MeV) = 
1.03 x 10 -5 ph cm~ 2 s _1 sr _1 and follows a power-law in energy with index 2.41 [19, 27]. 
This was determined by observing energies from 200 MeV to 100 GeV at angles above 10° in 
Galactic latitude (|6| > 10°). The DGRB flux accounts for roughly 25% of the total observed 
Fermi-LAT emission. 

Recently, Ref. [30] (hereafter FA12) measured the anisotropy of the DGRB. Taking the 
angular power spectrum of the DGRB in energies from 1 GeV to 50 GeV, they measure that 
the angular power above multipole moment I > 155 is nearly constant in I. This corresponds 
to angular scales of ~ 2° . The angular power observed in these high multipoles is significantly 
less than the noise due to finite measurement statistics expected for the Fermi-LAT. Therefore 
the Poisson term in the angular power spectrum is expected to come from unresolved sources 
rather than from measurement noise. For the DATA:CLEANED sample, the anisotropy is 
given in our Table 1. For many source classes, such as blazars, the DGRB anisotropy limits 
the total intensity of unresolved sources, which in turn constrains the contribution of that 
source class to the DGRB intensity. Here, we consider the contribution of unresolved blazars 
to the DGRB anisotropy and the corresponding contribution of such blazars to the DGRB 
intensity. 

Note that FA12 measured the anisotropy considering all the diffuse Fermi-LAT flux, 
including the Galactic diffuse which was removed in the DGRB measurement of Ref. [27]. 
Therefore, the intensity measured by FA12 differs from the DGRB spectral intensity in 
Ref. [27]. However, in stating the dimensionful anisotropy Cp, the Galactic diffuse emission 
is expected to only contribute at low £, so the value of Cp effectively neglects the Galactic 
diffuse emission and is consistent with the DGRB. Therefore, in this work, we compare 
calculations of the dimensionful Cp rather than the dimensionless Cp/(I) 2 . Also, because 

1 Note we use the term extragalactic diffuse gamma-ray background (DGRB) instead of the isotropic dif- 
fuse extragalactic gamma-ray background because the DGRB can have, and has been measured to have, 
anisotropies. 
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Table 1 



E (GeV) 


Blazar Cp 


DGRB C P 


1.04-1.99 
1.99-5.00 
5.00-10.4 
10.4-50.0 


4.52 x 10~ 18 
1.29 x 10~ 18 
8.89 x 10" 20 
2.36 x 10" 20 


4.62 x 10~ 18 
1.30 x 10~ 18 
8.45 x 10" 20 
2.11 x 10" 20 



Table 1. Comparison between the DGRB anisotropy Cp for our best-fit blazar model, given the 
constraints from the observed blazar source-count distribution and DGRB anisotropy, and, for com- 
parison, those observed in FA12. Values of Cp are given in (cm _2 s _1 sr _1 ) 2 sr. 



the anisotropy measurement does not depend on a model of Galactic diffuse emission, it has 
much smaller systematic errors than the DGRB intensity measurement. 

Here, we consider the consistency of blazars with the DGRB intensity and anisotropy. 
We begin by describing our general blazar model, an LDDE+SED-sequence model similar to 
that in ABH. By constraining this model with the blazar dN / dF and the DGRB anisotropy, 
we place limits on the LDDE model parameters. We place an upper limit on the contribution 
of blazars to the DGRB flux and predict the resolution of blazars by the Fermi-LAT at the 
end of 5 years of observations. We also place a lower limit on the blazar contribution to 
the DGRB anisotropy, which places a stringent constraint on all other sources of the DGRB 
intensity. Throughout the paper, we take a flat universe with the cosmological parameters 
tt m = 0.272, Q A = 0.728, and H = 70.2 km s" 1 Mpc" 1 [31]. Note, the use of h in the text 
refers to Planck's constant, and not the Hubble parameter. 

2 The LDDE+SED-Sequence Model of Blazar Intensity and Anisotropy 

2.1 The Blazar Gamma-ray Luminosity Function and SED-sequence 

The full model of blazars can be fully determined by two functions: the gamma-ray luminosity 
function (GLF) gives the comoving blazar density per unit luminosity and the SED gives the 
energy spectra of individual blazars. For a known redshift z and blazar luminosity L 7 (defined 
as vL v at hu = 100 MeV), the GLF is p^{L-y,z) and the SED is vL v (y, L<y). In this work, v 
denotes the blazar rest-frame frequency. 

By considering the blazar evolution as a function of both luminosity and redshift, we are 
effectively including both the FSRQ population, which has higher luminosity, and the BL Lac 
population, which has lower luminosity. Because FSRQs have higher luminosity, and therefore 
greater flux, than BL Lacs, the majority of unresolved blazars are BL Lacs. However, by 
considering all possible blazar luminosities and redshifts, we are including the contributions of 
both unresolved FSRQs and unresolved BL Lacs to the DGRB. The evolutionary properties 
of only the resolved FSRQs have been studied [32], though these blazars have been resolved 
and therefore do not contribute to the DGRB. 

Above a redshift of unity, the contribution of BL-Lacs to the blazar redshift distribution 
is negligible, so the two populations are expected to be the same. For redshifts greater than 
unity, our model successfully reproduces the redshift distribution of FSRQs in Ref. [33] , as 
is shown in figure 1. A Kolmogorov-Smirnov test (KS test) gives a probability of 78.6% that 
the distributions are consistent, indicating that the model reproduces the data well. 
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Figure 1. Shown is the best fit model line for all lFGL-observable blazars, both FSRQs and BL Lacs 
(solid line), relative to a measurement for the FSRQ population (dashed line) [33]. The consistency 
of the distributions is confirmed by a KS test, giving a probability of 78.6% that the distributions are 
consistent. Each distribution has been normalized to unity. . 

Because blazars and X-ray AGN are believed to be a related source population, observed 
at different projected angles, gamma-ray blazars should evolve similarly to X-ray AGN. For 
blazars, the gamma-ray jet is being observed, whereas for X-ray AGN, the X-ray emission 
from the disk is being observed. Hard X-ray AGN follow an X-ray luminosity function (XLF) 
which has been parameterized by Ref. [23]. The XLF model of Ref. [23] is a luminosity- 
dependent density evolution model, which has a peak evolution redshift which depends on 
luminosity, so AGN of different luminosities will evolve differently. 

The gamma-ray jet luminosity can be related to the X-ray disk luminosity Lx as 
P = 10 q Lx, where q is a scaling parameter and P = J L u dv is the bolometric blazar jet 
luminosity [26]. The bolometric luminosity from a blazar jet is proportional to the mass 
accretion rate m. For blazars with high accretion rates close to the Eddington limit, AGN 
X-ray disk luminosity also goes as Lx oc m [34-36]. Because black hole growth primarily 
takes place near the Eddington limit [37], the bolometric jet luminosity and the X-ray disk 
luminosity are proportional. If black hole growth were to take place far from the Eddington 
limit, the relative strengths of the bolometric jet luminosity and X-ray disk luminosity would 
vary [38]. 

In addition to the luminosity scaling parameter q, the GLF may have a different faint- 
end luminosity index 71 than the XLF. Also, the small opening angle of the blazar jet means 
that most AGN are not observed as blazars, so the fraction of AGN observed as blazars is 
parameterized as k. Note that this model does not require all AGNs to have blazar-like jets, 
but that the blazar population approximately follows that of the full AGN population, and is 
some fraction of it. In turn, a fraction of those can be seen as blazars. The comoving number 
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density per unit L 7 of gamma-ray blazars with luminosity L 7 at redshift z is then given as 

Pj(L 7 ,z) = K d ^p x (L x ,z) , (2.1) 

where px is the XLF. The precise relationship between X-ray AGN and gamma-rays blazars 
is not yet known. We are using the simple ansatz that they are related as shown in Eq. (2.1), 
as proposed by Inoue and Totani [26], and using physical data on gamma-ray blazars to 
constrain the three parameters. To the best of our knowledge, this model satisfies all current 
observations and constraints, and is therefore a viable possibility. 

The blazar SED sequence used in this work is that of Ref. [2-4] . This SED bins blazars 
by radio luminosity, giving vL v over nineteen orders of magnitude in energy for each of five 
radio luminosity bins. The blazar SED shows two peaks in the blazar luminosity: a peak in 
the radio associated with synchrotron emission and a peak in the gamma-ray associated with 
inverse Compton emission. A shift in these peaks with blazar luminosity indicates that the 
gamma-ray spectrum becomes harder for lower luminosity blazars. This has been observed 
by noting that FSRQs have high luminosity and soft spectra, while BL Lacs have harder 
spectra with lower luminosities [38-41]. The importance of this spectral hardening at low 
luminosities cannot be overstated: the gamma-ray spectra of unresolved, lower-luminosity 
blazars should be harder than the spectra of resolved blazars. 

To parametrize the SED, we use the model of Ref. [26], which fits the SED sequence of 
Ref. [2-4]. With this SED model, we can relate a blazar's bolometric luminosity to its gamma- 
ray luminosity to estimate the average energy spectrum for a blazar of that luminosity. The 
fully parameterized model of the SED may be found in Appendix 1 of ABH. 

2.2 LDDE Model Calculation of Blazar Number, Intensity, and Anisotropy 

For a blazar at redshift z with bolometric luminosity P, the observed gamma-ray flux in the 
energy band from E 1 ™ 111 to E™ 3 "* is 

S[z,P)= 2 / dv — exp 



hv 

-t | z, — — 
l + z 



(2.2) 



where di is the luminosity distance. The exp[— t(z, Eq)] factor accounts for photon absorption 
before reaching Earth. For our work, we choose the absorption factor of Gilmore et al. [42] , 
which has determined lower values of opacity than previous estimates and is consistent with 
the findings of Ref. [43]. For the special case of the blazar flux above 100 MeV (Eg 1111 = 
100 MeV and E^ ax = oo) we refer to the flux as F 100 . 

For a given GLF, the number count of blazars above a sensitivity F 7 is 

AT(>F 7 ) = 47r / dz— / dL 7 p 7 (L 7 ,z) , (2.3) 

JO dZ 7Z,lim( 2ii 7 7 ) 

where L 7 m is the gamma-ray luminosity below which a blazar at redshift z is no longer 
detectable for the sensitivity F 7 . We set the parameter z max = 5, but this does not affect the 
calculation significantly, since the GLF peak is at redshift of order unity. 

Integrating the flux of an individual blazar, Eq. (2.2), over all unresolved blazars, the 
total diffuse flux coming from unresolved blazars is 

dN _ 1 dz dX e -r( z ,E ) f L ^" z) dL ^JM Lu[E/h,P(Lj] ^ 



dEndAdtdQ. An J n dz Jt ' h E 

" u " ^.min 
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Here, E is the emitted photon energy, and Eq = E/(l + z) is the observed photon energy at 
Earth. A is area on Earth, t is time on Earth, and Q is solid angle in the sky. The integral 
over the comoving distance x accounts for the line-of-sight of the observation. Because the 
luminosity integral diverges if 71 > 1, we place a lower bound of L 7i min = 10 42 erg s^ 1 on 
this integral. This is an order of magnitude lower luminosity than any Fermi-LAT observed 
blazar [39, 40]. That is, we impose a step-function cutoff of blazar GLF at L 7i min- However, 
our best-fit value has 71 < 1, so the final result is not strongly dependent on this cutoff. 

The intensity per solid angle of the DGRB in the energy band from E™ 111 and E^ 1 ^ is 
given by 



The Poisson term of the blazar anisotropy per unit solid angle, Cp, is the square of the blazar 
flux, integrated over all unresolved blazars [25, 44]. For the blazar SED+GLF model, this is 
given by 



with V the comoving volume and S the single blazar flux in a particular energy band, from 
Eq. (2.2). This equation can be written as 



for luminosity distance di and comoving distance x- 

To calculate the threshold between observable and unobservable blazars, L 7 m , we use 
the index-dependent flux limit shown in figure 1 of FB10. This plot gives the flux limit as 
a function of blazar index for test statistic TS = 25, which is consistent with the DGRB 
flux [27] and anisotropy [30] analyses, which were done using the full Fermi-LAT 1-year 
catalog (1FGL). The blazar dN/dF of FB10 is tabulated for TS = 50 in the 1FGL, so for 
that calculation we use the flux limit shown in figure 1 of ABH. 

Recently, Ref. [45] has suggested another approach to determining the true source pop- 
ulation of blazars, to account for the observational biases discussed above. In that analysis, 
the correlation between the observed blazar flux and the observed blazar index distribution 
was estimated. By removing this observational bias, they determined the true blazar index 
distribution and dN/dF for the observed Fermi-LAT blazars. This method uses the index 
and flux of a given blazar to determine whether it is "resolved" or "unresolved," rather than 
the test statistic TS that has an intrinsically higher value for harder Fermi-LAT sources. 
In future analyses, such a technique may be used to limit the observational bias of the 
observed blazar source population. Using this method, Ref. [45] calculates an intensity of 
/(> 100 MeV) = 7.7+H x 10~ 6 ph cm 2 s 1 sr 1 for blazars, which is significantly larger 
than the value determined in this work when considering the DGRB anisotropy limits. In 
this work, we calculate the Fermi-LAT sensitivity threshold based on the TS values of FB10, 
as discussed above. 
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Figure 2. Shown is the best-fit model for the source-count distribution function dN/dF (solid line), 
given the constraints from the observed blazar source-count distribution and DGRB anisotropy. The 
data are from FB10. 

3 Results for Inferred Relation Between Blazar Intensity and Anisotropy 

Our adopted LDDE gamma-ray blazar model based on the X-ray AGN XLF requires three 
free parameters: q, 71, and k. All other parameters in the blazar model are fixed by other 
data, such as the AGN XLF and blazar SED sequence. 

First, we consider the consistency between the blazar dN/dF of FB10 with both the 
DGRB spectrum from Ref. [27] and the DGRB anisotropy from FA12. For no values of q, 
71, and k could the blazars account for both the full DGRB flux and DGRB anisotropy. In 
particular, values of the model which were consistent with the DGRB flux at high energies 
overproduced the DGRB anisotropy by a factor of ~ 5. The DGRB anisotropy constrains the 
blazar population much more than the DGRB intensity does. We also considered a model 
which was fit only to the blazar dN/dF. The best-fit model for this fit also overproduced 
the DGRB anisotropy, by a factor of 3.6. Therefore, the DGRB anisotropy itself places a 
significant constraint on the LDDE+SED-sequence blazar model. 

Considering only the consistency between the blazar dN/dF of FB10 with the DGRB 
anisotropy from FA12, we found that the DGRB anisotropy can be attributed entirely to 
unresolved blazars. The best-fit parameter values and errors to match the anisotropy and 
dN/dF are q = 4.2lig;Jjj, 71 = and log 10 (K/l(T 6 ) = 0.98±g;Jf . The model repro- 

duces well the blazar dN/dF and the DGRB C P , with a reduced x 2 /DOF = 0.76 for the 
combined data. To calculate the x 2 , we have assumed that each data point in the DGRB 
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anisotropy and each data point in the blazar dN/dF is independent, so we have included 
all 26 data points from the two data sets in our calculation. With an adopted Fermi-LAT 
5-year sensitivity of Fioo = 2 x 10~ 9 ph cm~ 2 s _1 [46], we expect 163ll34 5 blazars to be 
resolved, comprising 98.7 ± 0.1 percent of the total (resolved plus unresolved) blazar flux. 
The parameter values given by matching the dNf dF and anisotropy are significantly different 
from previous work which had found consistency between blazars and the DGRB flux (ABH 
and Ref. [26]), yet those models overproduce the DGRB anisotropy by a factor of ~ 5 (e.g. 
q = 4.19, 7i = 1.51, and log 10 (K/10~ 6 ) = 0.38 from ABH). 

The value of q indicates that the bolometric luminosity of a blazar jet is roughly ~16,000 
times more luminous than the X-ray from the accretion disk. The low value for 71 indicates 
that the anisotropy and flux due to blazars is dominated by high-luminosity blazars, as would 
be expected. The blazar fraction k indicates that of every ~105,000 AGN, one is observed as 
a blazar. If one assumes that every AGN is a blazar, this corresponds to an opening angle of 
~ 0.6 deg for the average blazar jet, similar to the intrinsic jet opening angle of ~ 1 deg [47]. 
Our fit to the DGRB anisotropy is given in Table 1. 

The dN/dF for our best-fit parameters is shown in figure 2. At low fluxes, the dN/dF 
flattens out instead of continuing upward to zero flux. This differs significantly from the 
broken power-law dN/dF assumed by much previous work. 

Given the blazar LDDE+SED-sequence model which best fits the blazar dN/dF and 
the DGRB anisotropy, we can infer the resulting constraints these measurements place on the 
flux of unresolved blazars. For the best-fit model, blazars only contribute 5.9^* % (68% CL) 
of the DGRB intensity above 1 GeV. Moreover, the blazar flux for our best-fit model does 
not produce greater than 8% of the DGRB flux in any single energy bin (68% CL). Above 
100 MeV, the best-fit model gives an intensity above 100 MeV of that is 4.0+^% (68% CL) 
of the DGRB. This is less than the value from Ref. [32], which predicts the intensity above 
100 MeV to be 9.3% of the DGRB from FSRQs alone. 

The spectrum of the inferred contribution of blazars to the DGRB relative to that 
observed is shown in figure 3. 

As discussed above, the best-fit blazar model is consistent with the full DGRB anisotropy. 
While it is clear that a valid blazar model cannot overproduce the DGRB anisotropy, it is 
possible that blazars are not a primary source of this anisotropy. To test this, we ran a 
further test of how little of the anisotropy is due to blazars. In this test, we included the 
three free parameters discussed above as well as a fourth free parameter which accounted 
for the fraction of the total measured Cp that is due to blazars. We found that in this 
test, the best-fit values of q, 71, and k were similar to the values discussed above for the 
combined dN/dF and DGRB anisotropy fit. However, we did find that the best-fit value 
of the blazar Cp was approximately the total DGRB Cp (98% of the total). This value of 
blazar Cp is consistent with the full DGRB anisotropy (68% CL), as is our best-fit model 
discussed above. In addition, this test finds that greater than 76% of the full DGRB Cp is 
due to blazars (68% CL). This provides a strong limit on any other potential contributor to 
the DGRB flux: any additional population which contributes to the remaining 92% of the 
DGRB flux must produce less than 24% of the DGRB anisotropy. 

The anisotropy-constrained LDDE+SED-sequence blazar model contributes much less 
of the intensity of the DGRB than in ABH. ABH had previously found consistency between 
the blazar dN/dF and the full DGRB intensity; here, blazars are constrained to less than 
5.9% of the DGRB intensity above 1 GeV. Though the model parameters q, 71, and k in 
ABH allowed for a blazar model which was consistent with the full DGRB intensity, that 
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Figure 3. The best- fit model for the total blazar intensity (solid red line), given the constraints 
from the observed blazar source-count distribution and DGRB anisotropy. The DGRB data are from 
Ref. [27]. 

model also severely overproduced the DGRB anisotropy. The DGRB anisotropy constraint 
changed the best-fit parameters significantly — e.g., the blazar faint-end index 71 changed 
from 1.51 (ABH) to 0.43 when the DGRB anisotropy was considered. This illustrates the 
importance of the DGRB anisotropy is constraining populations' contributions to the DGRB 
intensity. 

4 Discussion and Comparison to Previous Results 
4.1 Power-Law Extrapolation Models 

The blazar source-count distribution has been historically modeled by as a broken power-law 
in dN/dF when little was known regarding AGN population evolution and the background 
cosmological model. This simplifying assumption was adopted in FB10 and CKS. In such 
work, the broken power-law is fit to the source count distribution function above the point- 
source sensitivity observed by the Fermi-LAT. Below the point source resolution threshold, 
the faint-end power- law of this dN/dF continues to rise to zero flux, which produces the 
unphysical result of an infinite number of blazars within our observable horizon. To avoid 
this divergence, the dN/dF is expected to flatten at low flux, which occurs naturally in 
physically-derived models such as the LDDE+SED-sequence model. 

The spectra of individual blazars were approximated in FB10 as a distribution of single 
power-laws in energy. It was, however, shown in Refs. [39, 40] that a single power-law 
spectrum does a very poor job of fitting individual blazars, which display curvature and 
breaks in their spectra that are correlated with their total luminosity. Ref. [2-4] show that 
blazar gamma-ray spectra are soft for high-luminosity blazars and get much harder for low- 
luminosity blazars, and the wide variation of blazar spectra has been confirmed by Fermi- 
LAT [40]. The correlation of luminosity with spectrum leads to a hardening of the blazar 
spectra below threshold with respect to the spectrum above threshold. By not including a 



- 9 - 



flux-dependence in the blazar index distribution, FB10 has assumed a softer blazar spectrum 
below threshold than is expected from observations of low-luminosity blazars. 

Rather than considering the spectra of individual blazars, CKS considers only the ag- 
gregate blazar flux with a fixed power law in energy, which is inconsistent with the observed 
spectral variability of blazars and their dependence on luminosity [2-4, 40, 48]. This sim- 
plification can drastically affect the results, since the Fermi-LAT threshold flux significantly 
varies for blazars with spectral index, e.g, as in figure 1 of FB10. Therefore, a common flux 
threshold, as assumed in CKS, for all blazars independent of spectral index is an unnecessary 
and inaccurate approximation. The blazar intensity is most sensitive to unresolved blazars 
with hard spectra, including those with fluxes somewhat below threshold; conversely, the 
anisotropy is sensitive only to the blazar flux and is therefore most sensitive to blazars with 
flux just below and directly at the threshold. For an index-independent calculation of the 
flux threshold, the threshold used to calculate the intensity should be at a lower flux than 
that for the anisotropy calculation. Simple power-law blazar models such as those in CKS 
and FB10 have been historically used to extrapolate the blazar population below to fluxes 
below current observations. However, due to the reasons given above and elaborated on 
in Ref. [49], these extrapolation models can introduce significant unjustified approximations 
below the observational Fermi-LAT threshold and the behavior of the population near the 
threshold. 

4.2 Differences in Results Between the LDDE+SED Model and Power-Law 
Models 

The differences between the calculations using the LDDE+SED-sequence blazar model con- 
sidered here and the power-law models are threefold: the LDDE+SED-sequence gives a 
dN/dF which flattens below the Fermi-LAT threshold, the LDDE+SED-sequence model ac- 
counts for the hardening of the average blazar spectra for low-flux blazars, and the LDDE+SED- 
sequence model calculates an index-dependent Fermi-LAT threshold flux based on the harder 
low-flux blazars. Below Fi 00 ~ 4 x 10~ 9 ph cm" 2 s -1 , the blazar dN/dF for the LDDE+SED- 
sequence model rises less steeply than the broken power-law model of FB10. For comparison 
to CKS, we calculate dN/dS with the flux from 1-10 GeV, S, calculated as in Eq. (2.2) with 
£>min _ ^ Q e v and i?™ ax = 10 GeV. The comparison between these two models can be seen 
in figure 4. Although at high fluxes, where the Fermi-LAT data is well-constrained, the two 
models give similar source-count distributions, the LDDE+SED-sequence model flattens out 
at low flux instead of continuing to rise down to zero flux. This is expected in physically mo- 
tivated models, which have a continuous roll-off with flux rather than a hard break between 
two unrelated power-laws. 

The most striking feature of the LDDE+SED-sequence blazar model is due to the nature 
of the SED sequence itself. As blazar luminosity decreases, the peak of the inverse Compton 
emission shifts to higher energies. In the Fermi-LAT energy range, from 100 MeV to 100 
GeV, this trend is observable as a hardening in the intrinsic blazar spectrum at low fluxes. 

Where the Fermi-LAT point-source resolution is nearly complete, above -Fioo = 7 x 
10" 8 ph cm s , the LDDE+SED-sequence model matches the index distribution ob- 
served in FB10. This is shown in figure 5 by comparing the blazar index distribution of the 
LDDE+SED-sequence model above -Fioo = 7 x 10" 8 ph cm -2 s _1 to the measured blazar 
index distribution above the same threshold. A KS test gives a probability of 82.9% that the 
two index distributions above F\oo = 7 x 10~ 8 ph cm -2 s -1 are consistent. However, FB10 
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Figure 4. Shown is the best-fit model for the source-count distribution function dN/dS for flux S 
as in Eq. (2.2) with E™ ln = 1 GeV and _E™ ax = 10 GeV (solid red line), given the constraints from 
the observed blazar source-count distribution and DGRB anisotropy. For comparison, the best-fit 
dN/dS from FB10, used in CKS, is also shown (dashed black line). Note that at low fluxes, the two 
extrapolated models differ significantly. 

assumes that this index distribution is identical at all sensitivities while the LDDE+SED- 
sequence model accounts for the hardening in the blazar spectra at low flux. 

Additionally, there is a dependence of the threshold flux on the blazar index due to the 
bias of the Fermi-LAT towards observations of hard sources. For fluxes above 100 MeV, -Fioo; 
the flux threshold varies by two orders of magnitude with blazar index. For higher energy 
fluxes, this bias is much weaker, but still varies somewhat. For the flux from 1-10 GeV, for 
example, we calculate a threshold range from 2.4 — 7.4 x 10 _10 ph cm -2 s _1 for blazar indices 
from 1 to 1.5. These values are comparable to the value 3.6 x 10~ 10 ph cm -2 s _1 from FB10 
or the value 3.7 x 10 _10 ph cm~ 2 s~ x in CKS. 

For a blazar model with a single power-law spectrum, the blazar spectrum has a similar 
slope to the DGRB flux and therefore the intensity is a measure of the blazar flux in any 
bin. However, unresolved blazars have lower flux and therefore are expected to have a harder 
spectrum than the resolved blazars. For a source population with a harder spectrum than 
the DGRB, intensity is largely a measure only of the flux in the lowest energy bin, not a 
true indication of the source's contribution to the DGRB. Consider, for instance, the blazar 
model of ABH. The blazar population of ABH contributes only ~ 20% of the DGRB intensity 
because it has a much lower flux than the DGRB below 1 GeV. However, above 10 GeV, 
the ABH blazar model produces 100% of the DGRB flux. Therefore, DGRB flux, binned 
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Figure 5. Shown is the index distribution in our best-fit model above F 10 o — 7 x 1CP 8 ph cm~ 2 s _1 
(solid red line), given the constraints from the observed blazar source-count distribution and DGRB 
anisotropy. The measured index distribution above Fioo = 7 x 10~ 8 ph cm -2 s _1 from FB10 is shown 
for comparison (dashed black line). A KS test gives a probability of 82.9% that the two distributions 
are consistent. Both distributions in the figure have been normalized to unity. 

in energy, is a better determination of whether a particular model fits the DGRB than the 
integrated intensity is. 

Given the blazar LDDE+SED-sequence model which best fits the blazar dN/dF and 
the DGRB anisotropy, the LDDE+SED-sequence model constrains the flux from blazars to 
only contribute 5.9±? J% (68% CL) of the DGRB intensity above 1 GeV. This is a much 
stronger limit on the flux due to blazars than that from CKS or FA12, which derived limits 
of 30% and 19%, respectively. The blazar flux for the LDDE+SED-sequence best-fit model 
does not produce greater than 8% of the DGRB flux in any single energy bin (68% CL). 

5 Conclusions 

We have shown that the DGRB anisotropy can be composed entirely by blazar emission when 
adopting an LDDE+SED-sequence model. Given the constraints from the observed DGRB 
anisotropy, the blazars are likely only a sub-dominant contribution to the DGRB flux. That 
is, blazars cannot simultaneously account for the DGRB intensity as well as anisotropy in 
the LDDE+SED-sequence model we study. The LDDE+SED-sequence is an observationally- 
motivated physical-evolution model of the gamma-ray blazar population which successfully 
accounts for the shape of the blazar dN/dF below the current threshold, while taking into 
account the change in blazar spectra with flux, and the spectral dependence of the flux 
threshold. Models with a fixed source-count distribution power-law below a fixed Fermi-LAT 
point-source threshold fail to incorporate these crucial observed properties of the population. 

Using the LDDE+SED-sequence model, we find that unresolved blazars can account 
for all of the DGRB anisotropy, in general agreement with the conclusions of CKS and 
FA12. In the LDDE+SED-sequence model, blazars are expected to be the dominant source of 
anisotropy in the DGRB, comprising greater than 76% of the anisotropy (68% CL). However, 



- 12 - 



the constrained LDDE+SED-sequence blazar model produces only 5.7_ 1 ' % (68% CL) of the 
total DGRB intensity and less than 8% of the DGRB flux in any energy bin (68% CL). This 
is significantly less intensity than that calculated using the power-law models of FB10 or 
CKS. 

Other populations which contribute to the DGRB are highly constrained by these re- 
sults. Emission sources other than blazars, such as non-blazar AGN [50], millisecond pul- 
sars [51], star-forming galaxies [21, 52-55], radio galaxies [56], circum-galactic gas [57], or dark 
matter annihilation [44, 58-65], may be responsible for a significant portion of the DGRB 
flux, but such sources are constrained to provide less than 24% of the DGRB anisotropy. For 
many sources, this anisotropy constraint limits the amount of DGRB intensity those sources 
can provide. 

Recent studies have also shown that electromagnetic cascades from very high-energy 
emission of diffuse emission sources can contribute to the DGRB [66-69]. In some cases, this 
electromagnetic cascade can be used to constrain a particular source class's contribution to 
the DGRB intensity and anisotropy. In future analyses, such a technique may be used to 
further limit the blazar contribution to the DGRB in the LDDE+SED-sequence model. 

There are known deficiencies in the blazar LDDE+SED-sequence model's agreement 
with all currently observed blazar SEDs, and this must be addressed with further studies of 
the Fermi-LAT data. The blazar SED sequence considered here bases its gamma-ray spectral 
information on the few blazars observed by EGRET [2-4] . A new SED-dependence on blazar 
luminosity relation is strongly motivated, based on the latest Fermi-LAT results [39, 40, 48], 
and would likely greatly increase the accuracy of the predictions physical blazar population 
models. However, the blazar LDDE+SED-sequence is the best developed blazar-evolution 
model that goes far beyond the overly simple assumptions of a single power law for the full 
blazar population with a broken power-law source-count distribution. 

In summary, the LDDE+SED-sequence model of blazars is consistent with both the 
blazar source count distribution and the DGRB anisotropy. These results find that the 
blazar anisotropy limits blazars to be a sub-dominant contribution to the total DGRB flux. 
Physically-motivated models of the blazar SED-sequence and the blazar GLF will inevitably 
lead to a better understanding of blazars, non-blazar AGN, and the nature the contributors 
of the diffuse extragalactic gamma-ray sky. 

Acknowledgments 

JPH & KNA are partially supported by NSF CAREER Grant No. 11-59224. 
References 

[1] Fermi-LAT Collaboration Collaboration, M. Ackermann et. al., The Second Catalog of 
Active Galactic Nuclei Detected by the Fermi Large Area Telescope, arXiv: 1108. 1420. 

[2] G. Fossati et. al., Unifying models for X-ray selected and Radio selected BL Lac Objects, 
Astron. & Astrophys. 375 (Sept., 1997) 739-751, [astro-ph/9704113]. 

[3] G. Fossati et. al., A Unifying View of the Spectral Energy Distributions of Blazars, Mon. Not. 
Roy. Astron. Soc. 299 (1998) 433-448, [astro-ph/9804103]. 

[4] D. Donate et. at, Hard X-ray properties of blazars, Astron. & Astrophys. 375 (Sept., 2001) 
739-751, [astro-ph/0105203]. 



- 13 - 



[5] R. D. Blandford and A. Konigl, Relativistic jets as compact radio sources, Astrophys. J. 232 
(Aug., 1979) 34-48. 

[6] C. D. Dernier and N. Gehrels, Two Classes of Gamma- Ray -emitting Active Galactic Nuclei, 
Astrophys. J. 447 (July, 1995) 103-+. 

[7] P. Padovani et. al., Radio-loud AGN and the extragalactic gamma-ray background, Mon. Not. 
Roy. Astron. Soc. 260 (Feb., 1993) L21-L24. 

[8] F. W. Stecker, M. H. Salamon, and M. A. Malkan, The high-energy diffuse cosmic gamma-ray 
background radiation from blazars, Astrophys. J. Lett. 410 (June, 1993) L71-L74. 

[9] M. H. Salamon and F. W. Stecker, The blazar gamma-ray luminosity function and the diffuse 
extragalactic gamma-ray background, Astrophys. J. Lett. 430 (July, 1994) L21-L24. 

[10] J. Chiang et. al., The Evolution of Gamma-Ray-loud Active Galactic Nuclei, Astrophys. J. 452 
(Oct., 1995) 156-+. 

[11] F. W. Stecker and M. H. Salamon, The Gamma-Ray Background from Blazars: A New Look, 
Astrophys. J. 464 (June, 1996) 600-+, [9601120]. 

[12] J. Chiang and R. Mukherjee, The Luminosity Function of the EGRET Gamma-Ray Blazars, 
Astrophys. J. 496 (Mar., 1998) 752-+. 

[13] A. Miicke and M. Pohl, The contribution of unresolved radio-loud AGN to the extragalactic 
diffuse gamma-ray background, Mon. Not. Roy. Astron. Soc. 312 (Feb., 2000) 177-193. 

[14] P. Giommi et. al., Non-thermal Cosmic Backgrounds from Blazars: the contribution to the 
CMB, X-ray and gamma-ray Backgrounds, Astron. Astrophys. 445 (2006) 843-855, 
[astro-ph/0508034] . 

[15] T. Narumoto and T. Totani, Gamma-Ray Luminosity Function of Blazars and the Cosmic 
Gamma-Ray Background: Evidence for the Luminosity Dependent Density Evolution, 
Astrophys. J. 643 (2006) 81, [astro-ph/0602178]. 

[16] C. D. Dormer, Statistics of Cosmological Black Hole Jet Sources: Blazar Predictions for the 
Gamma-Ray Large Area Space Telescope, Astrophys. J. 659 (Apr., 2007) 958-975, [0605402]. 

[17] V. Pavlidou and T. M. Venters, The Spectral Shape of the Gamma-Ray Background from 
Blazars, Astrophys. J. 673 (Jan., 2008) 114-118, [arXiv: 0710. 0002]. 

[18] D. Bhattacharya, P. Sreekumar, and R. Mukherjee, Gamma-ray luminosity function of 
Gamma-ray bright AGNs, Res. Astron. Astrophys. 9 (2009) 85-94, [arXiv: 0811 .4388]. 

[19] Fermi-LAT Collaboration Collaboration, A. A. Abdo et. al., The Fermi-LAT high-latitude 
Survey: Source Count Distributions and the Origin of the Extragalactic Diffuse Background, 
Astrophys. J. 720 (2010) 435-453, [arXiv: 1003.0895]. 

[20] K. N. Abazajian, S. Blanchet, and J. Harding, The Contribution of Blazars to the Extragalactic 
Diffuse Gamma-ray Background and Their Future Spatial Resolution, Phys.Rev. D84 (2011) 
103007, [arXiv: 1012. 1247]. (ABH). 

[21] F. W. Stecker and T. M. Venters, Components of the Extragalactic Gamma Ray Background, 
Astrophys..]. 736 (2011) 40, [arXiv: 1012.3678]. 

[22] D. Malyshev and D. W. Hogg, Statistics of gamma-ray point sources below the Fermi detection 
limit, Astrophys. J. 738 (2011) 181, [arXiv: 1104. 0010]. 

[23] Y. Ueda, M. Akiyama, K. Ohta, and T. Miyaji, Cosmological Evolution of the Hard X-ray 
AGN Luminosity Function and the Origin of the Hard X-ray Background, Astrophys. J. 598 
(2003) 886-908, [astro-ph/0308140]. 

[24] G. Hasinger, T. Miyaji, and M. Schmidt, Luminosity-dependent evolution of soft X-ray selected 
AGN: New Chandra and XMM-Newton surveys, Astron. Astrophys. 441 (2005) 417-434, 



- 14 - 



[astro-ph/0506118]. 

[25] A. Cuoco, E. Komatsu, and J. Siegal-Gaskins, Joint anisotropy and source count constraints on 
the contribution of blazars to the diffuse gamma-ray background, arXiv: 1202.5309. 

[26] Y. Inoue and T. Totani, The Blazar Sequence and the Cosmic Gamma-Ray Background 
Radiation in the Fermi Era, Astrophys. J. 702 (2009) 523-536, [arXiv : 0810 . 3580]. 

[27] Fermi-LAT Collaboration Collaboration, A. A. Abdo et. ai, The Spectrum of the Isotropic 
Diffuse Gamma-Ray Emission Derived From First-Year Fermi Large Area Telescope Data, 
Phys. Rev. Lett. 104 (2010) 101101, [arXiv: 1002.3603]. 

[28] Y. Inoue et. al., The Cosmological Evolution of Blazars and the Extragalactic Gamma-Ray 
Background in the Fermi Era, arXiv : 1001 . 0103. 

[29] J.-Q. Xia, A. Cuoco, E. Branchini, M. Fornasa, and M. Viel, A cross- correlation study of the 
Fermi-LAT 7 -ray diffuse extragalactic signal, Mon.Not.Roy.Astron.Soc. 416 (2011) 2247-2264, 
[arXiv: 1103.4861]. 

[30] M. Ackermann, M. Ajello, A. Albert, L. Baldini, J. Ballet, et. ai, Anisotropics in the diffuse 
gamma-ray background measured by the Fermi LAT, arXiv : 1202 . 2856. 

[31] WMAP Collaboration Collaboration, E. Komatsu et. ai, Seven-Year Wilkinson Microwave 
Anisotropy Probe (WMAP) Observations: Cosmological Interpretation, Astrophys. J. Suppl. 192 
(2011) 18, [arXiv: 1001.4538]. 

[32] M. Ajello, M. Shaw, R. Romani, C. Dernier, L. Costamante, et. al., The Luminosity Function 
of Fermi- detected Flat-Spectrum Radio Quasars, Astrophys. J. 751 (2012) 108, 
[arXiv: 1110.3787]. 

[33] Fermi-LAT Collaboration Collaboration, A. A. Abdo et. al., The First Catalog of Active 
Galactic Nuclei Detected by the Fermi Large Area Telescope, Astrophys. J. 715 (2010) 429-457, 
[arXiv: 1002.0150]. 

[34] H. Falcke, E. Koerding, and S. Markoff, A Scheme to Unify Low-Power Accreting Black Holes - 
Jet- Dominated Accretion Flows and the Radio/X-Ray Correlation, Astron. Astrophys. 414 
(2004) 895-903, [astro-ph/0305335]. 

[35] A. Merloni, S. Heinz, and T. Di Matteo, A fundamental plane of black hole activity, Mon. Not. 
Roy. Astron. Soc. 345 (2003) 1057, [astro-ph/0305261]. 

[36] E. Gallo et. ai, A dark jet dominates the power output of the stellar black hole Cygnus X-l, 
Nature AUG (2005) 819-821, [astro-ph/0508228]. 

[37] A. Marconi et. al., Local Supermassive Black Holes. Relics of Active Galactic Nuclei and the 
X-ray Background, Mon. Not. Roy. Astron. Soc. 351 (2004) 169, [astro-ph/0311619]. 

[38] G. Ghisellini, L. Maraschi, and F. Tavecchio, The Fermi blazars' divide, Mon. Not. Roy. 
Astron. Soc. 396 (June, 2009) L105-L109, [arXiv : 0903 . 2043]. 

[39] Fermi-LAT Collaboration Collaboration, A. A. Abdo et. al., The Spectral Energy 

Distribution of Fermi bright blazars, Astrophys. J. 716 (2010) 30-70, [arXiv: 0912 . 2040]. 

[40] T. F.-L. Collaboration, Spectral Properties of Bright Fermi- detected Blazars in the Gamma-ray 
Band, Astrophys. J. 710 (2010) 1271-1285, [arXiv: 1001 .4097]. 

[41] B. Lott, E. Cavazzuti, S. Cutini, D. Gasparrini, and C. Dernier, General Properties of 
Fermi /LAT Active Galactic Nuclei, arXiv: 1205.3038. 

[42] R. C. Gilmore et. al., GeV Gamma-Ray Attenuation and the High-Redshift UV Background, 
Mon. Not. Roy. Astron. Soc. 399 (Nov., 2009) 1694-1708, [arXiv : 0905 . 1144]. 

[43] Fermi-LAT Collaboration Collaboration, A. A. Abdo et. at, Fermi Large Area Telescope 
Constraints on the Gamma-ray Opacity of the Universe, Astrophys. J. 723 (2010) 1082-1096, 



- 15 - 



[arXiv: 1005.0996]. 

[44] S. Ando, E. Komatsu, T. Narumoto, and T. Totani, Dark matter annihilation or unresolved 
astrophysical sources? Anisotropy probe of the origin of cosmic gamma-ray background, Phys. 
Rev. D75 (2007) 063519, [astro-ph/0612467]. 

[45] J. Singal, V. Petrosian, and M. Ajello, Flux and Photon Spectral Index Distributions of 
Fermi-LAT Blazars And Contribution To The Extragalactic Gamma-ray Background, 
Astrophysical Journal (2011) [arXiv: 1106.3111]. 

[46] LAT Collaboration, W. B. Atwood et. ai, The Large Area Telescope on the Fermi Gamma-ray 
Space Telescope Mission, Astrophys. J. 697 (2009) 1071-1102, [arXiv: 0902. 1089]. 

[47] A. Pushkarev, Y. Kovalev, M. Lister, and T. Savolainen, Jet opening angles and gamma-ray 
brightness of AGN, Astron. & Astrophys. 507 (Nov., 2009) L33-L36, [arXiv: 0910 . 1813]. 

[48] E. T. Meyer, G. Fossati, M. Georganopoulos, and M. L. Lister, From the Blazar Sequence to 
the Blazar Envelope: Revisiting the Relativistic Jet Dichotomy in Radio-loud AGN, 
Astrophys. J. 740 (2011) 98, [arXiv: 1107.5105]. 

[49] J. P. Harding and K. N. Abazajian, Comment on 'Joint Anisotropy and Source Count 
Constraints on the Contribution of Blazars to the Diffuse Gamma-Ray Background ', 
arXiv: 1204.3870. 

[50] Y. Inoue, T. Totani, and Y. Ueda, The Cosmic MeV Gamma-Ray Background and Hard X-Ray 
Spectra of Active Galactic Nuclei: Implications for the Origin of Hot AGN Coronae, Astrophys. 
J. Lett. 672 (Jan., 2008) L5-L8, [arXiv: 0709. 3877]. 

[51] C. A. Faucher-Giguere and A. Loeb, The Pulsar Contribution to the Gamma-Ray Background, 
JCAP 1001 (2010) 005, [arXiv: 0904. 3102]. 

[52] B. D. Fields, V. Pavlidou, and T. Prodanovic, Cosmic Gamma-Ray Background from 
Star-Forming Galaxies, Astrophys. J. 722 (2010) L199, [arXiv: 1003.3647]. 

[53] N. Chakraborty and B. D. Fields, Inverse Compton Contribution to the Star-Forming 
Extragalactic Gamma- Ray Background, arXiv: 1206.0770. 

[54] B. C. Lacki, S. Horiuchi, and J. F. Beacom, The Star-Forming Galaxy Contribution to the 
Cosmic MeV and GeV Gamma-Ray Background, arXiv: 1206.0772. 

[55] Fermi LAT Collaboration Collaboration, M. Ackermann et. ai, GeV Observations of 
Star-forming Galaxies with Fermi LAT, Astrophys. J. 755 (2012) 164, [arXiv: 1206. 1346]. 

[56] Y. Inoue, Contribution of the Gamma-ray Loud Radio Galaxies Core Emissions to the Cosmic 
MeV and GeV Gamma-Ray Background Radiation, Astrophys. J. 733 (2011) 66, 
[arXiv: 1103.3946]. 

[57] R. Feldmann, D. Hooper, and N. Y. Gnedin, Circum- Galactic Gas and the Isotropic Gamma 
Ray Background, arXiv : 1205 . 0249. 

[58] S. Ando and E. Komatsu, Anisotropy of the cosmic gamma-ray background from dark matter 
annihilation, Phys. Rev. D73 (2006) 023521, [astro-ph/0512217]. 

[59] F. Miniati, S. M. Koushiappas, and T. Di Matteo, Angular Anisotropics in the Cosmic 
Gamma-ray Background as a Probe of its Origin, Astrophys. J. 667 (2007) L1-L14, 
[astro-ph/0702083] . 

[60] E. A. Baltz et. ai, Pre-launch estimates for GLAST sensitivity to Dark Matter annihilation 
signals, JCAP 0807 (2008) 013, [arXiv: 0806. 2911]. 

[61] J. M. Siegal-Gaskins, Revealing dark matter substructure with anisotropics in the diffuse 
gamma-ray background, JCAP 0810 (2008) 040, [arXiv : 0807 . 1328]. 

[62] J. M. Siegal-Gaskins and V. Pavlidou, Robust identification of isotropic diffuse gamma rays 



- 16 - 



from Galactic dark matter, Phys. Rev. Lett. 102 (2009) 241301, [arXiv: 0901 . 3776]. 

[63] B. S. Hensley, J. M. Siegal-Gaskins, and V. Pavlidou, The detectability of dark matter 

annihilation with Fermi using the anisotropy energy spectrum of the gamma-ray background, 
Astrophys. J. 723 (2010) 277-284, [arXiv: 0912. 1854]. 

[64] M. Fornasa, L. Pieri, G. Bertone, and E. Branchini, Anisotropy probe of galactic and 

extra-galactic Dark Matter annihilations, Phys. Rev. D80 (2009) 023518, [arXiv : 0901 . 2921]. 

[65] K. N. Abazajian, S. Blanchet, and J. Harding, Current and Future Constraints on Dark Matter 
from Prompt and Inverse-Compton Photon Emission in the Isotropic Diffuse Gamma-ray 
Background, Phys. Rev. D85 (2012) 043509, [arXiv: 1011 . 5090]. 

[66] T. M. Kneiske and K. Mannheim, BL Lac Contribution to the Extragalactic Gamma-Ray 
Background, Astron. Astrophys. 479 (Feb., 2007) 41-47, [arXiv: 0705. 3778]. 

[67] T. M. Venters, Contribution to the Extragalactic Gamma-ray Background from the Cascades of 
Very-high Energy Gamma Rays from Blazars, Astrophys. J. 710 (2010) 1530-1540, 
[arXiv: 1001. 1363]. 

[68] T. M. Venters and V. Pavlidou, Probing the Intergalactic Magnetic Field with the Anisotropy of 
the Extragalactic Gamma-ray Background, arXiv: 1201.4405. 

[69] Y. Inoue and K. Ioka, Upper Limit on the Cosmological Gamma-ray Background, Phys. Rev. 
D86 (2012) 023003, [arXiv: 1206.2923]. 



-17- 



